Dynamics of two trapped Brownian particles: 
shear-induced cross-correlations 



O 
(N 

Oh: 



Jochen Bammert, Lukas Holzer, Walter Zimmermann 
Theoretische Physik I, Universitdt Bayreuth, D-95440 Bayreuth 
(Dated: Received: September 8, 2010/ Revised version: September 8, 2010) 

The dynamics of two Brownian particles trapped by two neighboring harmonic potentials in a 
linear shear flow is investigated. The positional correlation functions in this system are calculated 
analytically and analyzed as a function of the shear rate and the trap distance. Shear-induced 
cross-correlations between particle fluctuations along orthogonal directions in the shear plane are 
found. They are linear in the shear rate, asymmetric in time, and occur for one particle as well as 
between both particles. Moreover, the shear rate enters as a quadratic correction to the well-known 
correlations of random displacements along parallel spatial directions. The correlation functions 
depend on the orientation of the connection vector between the potential minima with respect 
to the flow direction. As a consequence, the inter-particle cross-correlations between orthogonal 
fluctuations can have zero, one or two local extrema as a function of time. Possible experiments for 
detecting these predicted correlations are described. 

PACS numbers: 05.40.-a,33.15.Vb,47.15.G- 
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I. INTRODUCTION 

The Brownian dynamics of particles in fluids is of high 
relevance in many fields of natural and applied sciences. 
It is strongly affected by the interplay between the par- 
ticles via the liquid, the so-called hydrodynamic interac- 
tion. Especially in the field of microfluidics this nonlin- 
ear interaction plays an important role in subjects such 
as Taylor dispersion [l| or fluid mixing 0, Q . In a quies- 
cent fluid there is already a considerable understanding 
of the dynamics of Brownian particles and their interac- 
tions [3, Q . Investigations on the positional correlation 
functions of two trapped particles give further insight in 
the coupling between thermal motion and the hydrody- 
namic interaction among them @, 01- However, the un- 
derstanding of this interplay in typical laminar flows, like 
in a linear shear flow or in a Poiseuille flow, is still far 
from complete although it is the origin of a number of 
interesting phenomena. For example, polymers exhibit 
in shear flows already at small values of the Ileynolds 
number the so-called molecular individualism d, Q, elas- 
tic turbulence, and spectacular mixing properties in the 
dilute regime [lO[ . The complex interplay between Brow- 
nian motion and hydrodynamic interaction also affects 
considerably the conformational distribution functions of 
tethered polymers in flows and their dynamics [ll| - [l9| . 
When polymers are attached to a wall and subjected to 
a flow, an additional time-periodic behavior influences 
the dynamics [20l - |2^ . which shows similar features as 
tumbling polymers in shear flows [1, [25l - [27| . 

Recent theoretical investigations on the fluid-velocity 
fluctuations in shear flows show that in contrast to 
quiescent fluids or uniform flows, cross-correlations be- 
tween velocity fluctuations along and perpendicular to 
the streamlines occur [2^ [2^ . For free Brownian par- 
ticles in a linear shear flow in x direction, where the 
shear plane is parallel to the xy plane, one also expects a 
cross-correlation between the orthogonal positional fluc- 



tuations X and y of the particles, i.e. {xy) ^ [30l - l32| . In 
this case, random jumps of a particle between neighbor- 
ing parallel streamlines lead to a change of the particle's 
velocity. For example, a positional fluctuation, y, per- 
pendicular to the streamlines may cause a fluctuation, 
i, along the streamlines and contributes in this way to 
the correlation function (iy), which reflects the shear- 
induced coupling between fluctuations along orthogonal 
directions. 

The theoretical considerations described in Ref. [s^ 
show that shear-induced cross-correlations between per- 
pendicular random particle displacements, like {xy) ^ 0, 
survive if a particle experiences some constraints such 
as a harmonic potential. This is important from various 
points of view. First, these cross-correlations are inher- 
ently present in bead-spring models, which are used to 
describe polymer dynamics in shear flows, because the in- 
dividual beads along the chain are bound to their neigh- 
bors. Second, this knowledge facilitates the experimen- 
tal detection of these cross-correlations, because mea- 
surements of particle fluctuations in the spatially limited 
area of the trapping potential can be performed in a con- 
trolled manner compared to tracing free Brownian par- 
ticles. According to this strategy, the cross-correlations 
{xy) of trapped particles in a linear shear flow have been 
measured directly for the flrst time and the results are in 
good agreement with the theoretical predictions [13] ■ 

The optical tweezer technique, employed in Ref. [s^ ]. 
triggered a number of further direct observations of par- 
ticle fluctuations. These include inspiring studies on sin- 
gle polymers pll- [la ^35|. the propagation of hydrody- 
namic interactions [36| . wall effects on Brownian mo- 
tion [13, HI], t wo-p oint microrheology [1^, particle sort- 
ing techniques j40l - |43j , the determination of the effective 
pair potential in colloidal suspensions [44 k and many 
other investigations in microfluidics, cf. 4^ |4^. By 
femto-Newton measurements anti-correlations have been 
detected between two hydrodynamically interacting and 
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neighboring particles in a quiescent fluid, each one cap- 
tured by a laser-tweezer potential 0. This was recently 
extended in Ref. [s^], where shear- induced inter-particle 
anti-correlations between orthogonal motions of the two 
particles have been found. 

The present work focuses on the question, which kind 
of cross-correlations can be expected between hydrody- 
namically interacting particles in linear shear flows. Such 
inter-particle correlations along a single polymer and be- 
tween different polymers in flows influence their dynam- 
ics. This paper is an extension of the theoretical work on 
the single particle dynamics described in Ref. to a 
pair of two hydrodynamically interacting point-particles 
with an effective hydrodynamic radius and trapped in a 
linear shear flow. It provides the theoretical background 
for the experimental results on the two-particle correla- 
tions presented in Ref. [s^l • The correlation functions be- 
tween the different particle displacements arc calculated 
analytically and we show how a second Brownian particle 
influences the stochastic motion and the positional prob- 
ability distribution of its neighbor compared to the single 
particle case [H, [131 • In addition, we find that the anti 
cross-correlations between two fluctuating particles in a 
quiescent fluid, as described in Ref. 0, experience shear- 
induced corrections. We also describe the occurrence of 
shear-induced anti cross-correlations between the fluctua- 
tions of the two particles along two orthogonal directions. 
The results depend significantly on the orientation of the 
connection vector between the two traps with respect to 
the fiow direction. We focus on the leading order con- 
tributions to the correlation functions and neglect the 
effects of finite size and rotations of the particles. 

The structure of the paper is as follows: In sect, [ill the 
model equations are introduced and their formal analyti- 
cal solution is presented, which results in the calculation 
of the correlation functions for the random particle dis- 
placements. In sect, mil the results for three representa- 
tive setups of the two-particle system are discussed in de- 
tail, where the connection vector between the two poten- 
tial minima is cither parallel, perpendicular, or oblique to 
the shear-fiow direction. In addition, the results are com- 
pared with direct simulations of the Langevin equation 
for representative examples. For the parallel few 
experimental and theoretical results have already been 
described in Ref. [131, where a good agreement between 
experiment and theory was found. The article closes with 
a discussion and further possible applications in sect. II VI 



II. EQUATIONS OF MOTION AND THEIR 
SOLUTION 




b) ^ = 7r/2 
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FIG. 1. Two Brownian particles are kept by linear spring- 
forces iY,2 close to the minima of two corresponding harmonic 
potentials at qi,2. Both particles are simultaneously exposed 
to a linear shear flow u(r). The angle between the flow 
direction u(r) and the vector qi2 = qi — q2 connecting the 
potential minima at the distance d — |qi2| is either zero, as 
in part a), or (j) = n/2 as in b), or </!> = n/4 as in c). 



are discussed briefly. 



A. Model equations 

We consider two Brownian point-particles with effec- 
tive hydrodynamic radius a, immersed in a Newtonian 
fluid of viscosity 77 at the positions r.^ = {xi,yi,Zi) with 
i = 1,2. Each particle is held by a linear restoring force. 



(1) 



close to the minimum of a corresponding harmonic 
potential. 



2 - q*) 



(2) 



with the spring constant k. The distance between the po- 
tential minima is labeled with d. Uncharged polystyrene 
latex beads of micrometer size that are trapped by laser 
tweezers experience such a potential as described by 
eq. ^ %m^. 

Both trapped particles are exposed to a linear shear 
flow in the x direction with the shear plane parallel to 
the xy plane and the shear rate 7: 



The basis of our investigation is a Langevin model 
that describes the over-damped dynamics of two parti- 
cles, each held by a harmonic potential in a linear shear 
flow. In this section the equation of motion is introduced 
and solved in order to calculate the correlation functions 
analytically. They consist of different eigcnmodes, which 



u(r) = 7ye^ , 



(3) 



This flow causes a drag force on the hydrodynamically 
interacting Brownian particles, which is in competition 
with the restoring force (H)). In a recent experiment the 
interplay between these two forces has been studied [3J] , 
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where the size of the particles was about 5/im and the 
shear rate about 50s~^. At the length scale of the ex- 
cursions of the particles from their potential minima and 
the distance between the two beads, the Reynolds num- 
ber is small and therefore, we can describe the fluid mo- 
tion around the beads in the Stokes limit. Consequently, 
the over-damped particle dynamics is described by the 
Langevin equations. 



V I 



V, = n{v,) +'Hjjfj 
with the four 3x3 mobility matrices 
1 



c 



n.12 — rt.21 — -:- 



C4ri2 



ri2 ® ri2 



(4) 

(5a) 
(5b) 



including the Stokes friction coefficient Q = Qirrja for a 
single point-particle. H^. The latter two matrices de- 
scribe the hydrodynamic interaction between two point- 
particles in terms of the Oseen tensor Q and X represents 
the unity matrix. The dyadic product (tensor product) ® 
has been used as well as the distance vector ri2 := ri — r2 
between the beads, with ri2 = |ri2|- 

The stochastic forces in the Langevin model have their 
origin in the velocity fluctuations of the surrounding liq- 
uid. In a quiescent fluid these random forces are un- 
correlated along orthogonal directions in the bulk jsof . 
This assumption is kept in our Langevin model, because 
shear-induced cross-correlations between the stochastic 
forces along orthogonal directions are expected to be 
small [UlllilEl]. So, for the contribution ff (t) in cq. ^ 
we assume a zero mean and a vanishing correlation time 

i: 

(ff(i))=0, (6a) 
(ff (i) ® f/(t')) = 2kBTnrj6{t ~ t') . (6b) 

The strength of the stochastic forces is proportional to 
the thermal energy ksT. 

The orientation of the connection vector between the 
potential minima, qi2 := qi — q2, with respect to the flow 
direction is described by the angle (j). It has a strong in- 
fluence on the correlation functions of the positional fluc- 
tuations of the particles. For this reason, we investigate 
three characteristic setups, where qi2 is either parallel to 
the external flow u (<^ = 0), or perpendicular (cj) = 7r/2), 
or oblique (0 = 7r/4) as sketched in fig. [1] 



B. Solutions and relaxation times 

There are two characteristic time scales in the system. 
One is determined by the inverse shear rate 7"^ and the 
other one is given by the relaxation time r := (/k of the 
particles in the two identical potentials. Their ratio gives 
the dimensionless Weissenberg number. 



Wi := jr , 



(7) 



which will be useful for the further discussion. 

The first step in the solution of cq. ([4]) is to rewrite 
the equation of motion in a more compact form by intro- 
ducing the positional vector R = (ri,r2) with six com- 
ponents and the 6x6 mobility matrix "H, 



n 



T~Li2 T~L 
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(8) 



composed of the sub-matrices 'Hij from eqs. ([5]). The 
shear flow in eq. ^ can be written in an analogous man- 
ner with the 6x6 shear rate tensor U, 



U(R) =Z^R, 

with hli2 = W45 = 7 and all other lAki 
of motion (|4]) then takes the form, 

R = WR + fc-H (Q - R) 



(9) 

0. The equation 
F, (10) 



where the vector Q — (qi,q2) describes the positions of 
the two potential minima being separated by the distance 
d. The stochastic contribution F in eq. (|10|) is obtained 
from eqs. ([S]): 



(F(t) 



(F(t)) 
F(t')) 



0, 

2kBTnS{t 



t'). 



(11a) 
(lib) 



We assume two well separated point-particles with small 
values of a/d and small fluctuations around their mean 
positions, i.e. ksT / ika^) ^ 1. The experimental results 
described in Ref. [SJ], which were obtained for a/d ~ 1/4 
and a magnitude of the fluctuations below a/10, are well 
described within this approximation. 

Since wc investigate the particle fluctuations, the mean 
position R^ := (R(i)) has to be determined flrst. In the 
case (j) = the mean positions are identical with the 
locations of the potential minima: Rg = Q. For cj) = it /2 
and = 7r/4, R^ is obtained numerically by determining 
the stationary solution of eq. pO|) in the absence of noise. 
Disregarding the hydrodynamic interaction between the 
two particles one flnds the analytical expressions. 



Rw2= o(Wi,l,0,-Wi,-l,0) , 



(12a) 



R^4 = - (\/2 + Wi, V2, 0, - V2 - Wi, -72, 0) , 



(12b) 

which may serve as an approximation of the station- 
ary solution and as the starting point of the numerical 
iteration. 

The equation of motion for the particle fluctuations 
R = (xi, yi, ^1, 2^2, 1/2, ^2) arc obtained by the ansatz R = 
R^ -I- R and the linearization of eq. (fTU)) with respect to 
R: 

R = Z^R - fc-HR + /c [V-HR] (Q - R0) + F. (13) 
Here the mobility matrix T-L and its derivative are eval- 



uated at the exact mean positions Ra. For 
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one has = Q and the third contribution on the 
right hand side vanishes. Introducing the matrix K, := 
[V (8) ('H(Q — R^))]^, eq. (fT3|) can be rewritten to 



a) 



V2 



b) 



R=_(fc-H-W-fc/C)R + F = -XR + F, (14) 
and this linear equation has the formal solution: 

R(t) = e-*-^R(0) + / dt'e(*'-*)-^F. (15) 
Jo 



By introducing the scaled deviation 

~ R . , „ IkeT 

R = — ■= , with B — , 

x/B k 



(16) 



and taking into account the statistical properties of the 
stochastic forces as given by eq. pT|) . one can determine 
by a straight-forward calculation, assisted by computer 
algebra, the correlation matrix C{t) defined by, 



C(t) {R{0)(g)'k{t)) for t>0. 



(17) 



The brackets (•) denote the ensemble average over a large 
number of particle trajectories. The elements Cki{t) of 
the Matrix C{t) can be represented as a sum of six expo- 
nentially decaying contributions. 



(a, k, I = 1...6) , 



(18) 



where the coefficients ga.ki depend on the Weissenberg 
number Wi ~ 7T and on the distance between the po- 
tential minima d. 

The origin of the relaxation times l/i?e(AQ), given by 
the eigenvalues of the matrix A4 := kH — U — kIC, 
can be explained as follows: After a stochastic kick that 
pushes the particles away from their mean positions, the 
potential forces start to pull them back. During this re- 
laxation the particle motion can be decomposed into par- 
allel or anti-parallel translations as illustrated in fig. [21 
It is the hydrodynamic interaction between the particles, 
which accelerates or damps this process, since the result- 
ing hydrodynamic forces depend on the relative particle 
motions. The beads relax faster, if they move in the same 
direction. The whole relaxation process is described by 
six relaxation rates, Aq, two for each spatial direction 
and in the two cases = 7r/2 and = 7r/4 some of them 
may even be complex. For the configuration = 0, cf. 
fig.[T]a), two relaxation modes coincide, so there are only 
four instead of six different relaxation times. 



III. CORRELATION FUNCTIONS 

In the two-particle system, the fluctuation statistics 
of one particle is influenced by its neighbor. We call 



FIG. 2. The motion of particle 1 causes via the hydrody- 
namic interaction a force f^' on particle 2 and vice versa. 
In part a), where the two particles move in the same direc- 
tion, the resulting hydrodynamic forces accelerate the motion. 
This situation corresponds to the parallel relaxation. Part b) 
shows the anti-parallel case, where the forces fj^' decelerate 
the motion. 



the corresponding correlation functions 'one-particle cor- 
relations' and for the inter-particle cross-correlations be- 
tween the random displacements of two different particles 
we use the notation 'inter-particle correlations'. 

In a previous study on hydrodynamic interactions be- 
tween two trapped Brownian particles in a quiescent 
fluid, anti cross-correlations between their random mo- 
tions along parallel spatial directions were found Q . We 
show, that the shear flow alters these correlations and ad- 
ditionally induces inter-particle cross-correlations along 
orthogonal directions in the shear plane. 

In this section the exact expressions and the approx- 
imations of the correlation functions Cki{t), as defined 
by eq. pT]) . are discussed in detail and they are also 
compared with numerically obtained solutions of the 
Langevin equation (fTO|) . The behavior of Cki{t) depends 
on the trap distance d, and in the limit d — > 00, our 
formulas become identical to the recently presented re- 
sults for a single trapped particle in a linear shear flow 
[33I ) . Since our results depend on the angle between the 
connection vector qi2 and the flow direction the three 
characteristic configurations, as sketched in fig. (TJ are 
analyzed. 



A. Parallel case: = 

At first, we consider the two-particle configuration 
with the connection vector qi2 parallel to the flow lines 
u as sketched in fig.[T]a), e.g. Q = d/2(l, 0, 0, -1, 0, 0) . 
The discussion of the one-particle correlations in sect. 
nil A II is complemented by the analysis of the inter- 
particle correlations in sect. IIII A2I 

Similar to the case of two trapped particles in a qui- 
escent fiuid in Ref. Q there are four different relaxation 
rates describing the four relaxation times in the system, 
cf. secLUm 



Ai = 



A, = 



1 + 2^ 

T 



A. = 



A4 = 



1 - 2/i 

7 

r 



(19a) 
(19b) 



The parameter < := 3a/ (4(i) < 3/8 is a measure for 
the distance between the traps. Ai and A3 correspond to 
the particle motions parallel and anti-parallel to the con- 
nection vector qi2 (longitudinal displacements), while A2 
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FIG. 3. The autocorrelations, {xi{0)xi{t)) = {x2{0)S:2{t)), 
along the flow direction are shown for Wi = (dash-dotted 
line) and for Wi = 1/2 (dashed line). The correlation 
functions perpendicular to the flow direction, (yi{0)yi{t)) ~ 
{zi{0)zi{t)) with j = 1,2 (solid line), do not depend on the 
Weissenberg number. All curves are obtained for a distance 
d — Aa and t is given in units of r. 



and A4 belong to the particle relaxations perpendicular 
to qi2 (transversal displacements). 



1. One-particle correlations 

The autocorrelations are identical for both particles 
but they are different for the longitudinal and transver- 
sal displacements: {xi{0)xi{t)) = {x2{0)x2{t)) and 

(yi(0)yiW> = (^2 (0)^2(0) = {h{0)h{t)) = (f2(0)z2W>. 

The expressions 



(ii(0)ii(t)) = i(e-^^* + e-^^*) 



Wi f-{l + ^i)e-^' 



Wi^ / (1 - n)e 



(yi(O)yi(i)) = (zi(O)zi(t)) ^ \ {e'^^' + e"^^*) 



(20a) 
(20b) 



are exponentially decaying in time and both functions 
are plotted in fig. [3] for different values of the Weis- 
senberg number Wi. Due to the scaling (fT5|) the value 
of (yi(O)yi(O)) is 1/2. The correlation functions de- 
pend via /i on the trap distance d , which leads to inter- 
esting corrections to the autocorrelations compared to 
the case of one isolated particle in Ref . [s^ . 

The distinct relaxation rates given by eqs. (fT9)) lead to 
different autocorrelations of particle displacements along 
and perpendicular to qi2, independent of the parame- 
ter Wi. This is also indicated in fig. [3] by the difference 



FIG. 4. Shear-induced cross-correlations of a single particle 
in the parallel case for Wi = 1/2 and d = Aa respectively 
d = lOOo. {x\{t)y\{0)) has a maximum around t' « r, which 
does not depend on Wi and only very weakly on d. 



between the correlations {xi[Q)xi{t)) (dash-dotted line) 
and (j/i(0)yi(t)) (solid line). The latter one is indepen- 
dent of the Weissenberg number, which is similar to the 
case of an isolated trapped particle, where only the auto- 
correlation in fiow direction (|20ap depends on Wi^ [s^l . 

Any translation of a particle in y direction is coupled 
via the flow profile ([3]) to a change of the particle's veloc- 
ity in X direction, which leads to a change of the parti- 
cle's positional fluctuation along the x direction. Conse- 
quently, the particle fluctuations in the x and y directions 
become correlated due to the shear flow. The result- 
ing cross-correlations between the displacements along 
orthogonal directions in the shear plane are linear func- 
tions of the par ameter Wi, similar to the single particle 
case in Refs. ]33l[34j . However, compared to these results, 
we obtain for the two-particle system an additional de- 
pendence on the trap distance d. The time-dependence 
of the cross-correlations is given by the following expres- 
sions: 



-Aot 



(yi(0)£i(t)) = ^(e- 



-Aot 



2 - 3/x 



Wi / (I - /f)e-^^* (l-F/i)e-^^* 
2^ V 2-3^ 2-f 3^ 



(21a) 



(21b) 



Both functions are plotted in flg. [4] for two different dis- 
tances. The two sets of curves indicate that the mag- 
nitude of the cross-correlation increases weakly with de- 
creasing values of d. The cross-correlation functions in- 
volving the z coordinate vanish as in Ref. [33| . 

The time-asymmetry of the shear-induced cross-corre- 
lations, namely {yi{Q)xi{t)) 7^ (xi(O)yi(t)) with t > 0, 
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can be explained in the following way: A random particle 
displacement at t = in y direction leads immediately 
after the kick to a linear growth of the particle's x coordi- 
nate due to the larger flow velocity u(r) at a larger value 
of the y coordinate. Consequently, for small values of t, 
the product yi{0)xi{t) grows in time until the particle 
is pulled back by the linear spring force, which happens 
on the time scale t = ^//c. The result is a maximum 
in the correlation function (j/i(0)ii(t)) at a time of the 
order of r. Considering the effect of a random kick in x 
direction at t = 0, the particle does not jump between 
streamlines of different velocity and therefore the cross- 
correlation (ii (0)yi (t)) does not show this maximum and 
decays exponentially. 

By replacing the time t by ~t in cq. ()21ap the two func- 
tions (|21a|) and (j21bp can be combined to one correlation 
hmction {xi(0)yi{t)) , where t can now take positive and 
negative values. This function is asymmetric with re- 
spect to time-reflections t — —t. A similar behavior was 
previously found for the fluctuations of the fluid- velocity 
in a shear flow j28| . 

The static correlation functions given by eqs. pOa|l - 
(j21bp for i = determine also the positional distribution 
function P(r) of a Brownian particle in a potential as 
described in more detail in Ref. [l^. It is an interesting 
question, how the single particle distribution in a shear 
flow is changed by the presence of a second one. 

In a linear shear flow P{r) has an elliptical shape in 
the shear plane and the angle 9 between the major axis 
of the ellipse and the flow direction is determined by the 
equation 



tan 6* = - 
2 



[xiyi) 



with 



G = 



\{S:iyi)\ 



V4 + G2 -G 



{xiyi) 



(22) 
(23) 



Using the d-dependent static correlations from eqs. ([20 
and (1211) we obtain 



G = Wi 



. l + 3/i2 



1 - V 



(24) 



Consequently the inclination angle is a function of the 
parameter d too. In fig.[5]tan(0) is shown as a function of 
jjL = ia/ (Ad) for three different values of the Weissenberg 
number Wi and in all three cases tan(0) decreases consid- 
erably when the two particles approach. Since 6 changes 
in the same manner by decreasing d or increasing Wi, 
the shear-flow effects can be considered to be amplified 
by the presence of the second particle. 



2. Inter-particle correlations 

The motion of the two trapped Brownian particles is 
coupled via the hydrodynamic interaction. In a quiescent 
fluid this coupling leads to a cross-correlation between 




FIG. 5. The inclination angle 9 of the single particle distri- 
bution is shown as function of /x = 3a/(4d) for three different 
values of the Weissenberg number Wi. 



the thermal fluctuations of the two particles along the 
same direction Q. Since this cross-correlation is nega- 
tive as a function of time, cf. fig. [6l the notion 'anti 
cross-correlation' is used. For the anti cross-correlation of 
the longitudinal displacements we obtain a shear-induced 
correction, similar as for the one-particle autocorrelation 
in cq. ()20ap . which is proportional to Wi^: 



(xi(0)£2(i)) ^-Ae 



-{l + ^l)e 



-Alt 



-Alt 



-As* 



Aat 



Wi^ /-(I - ^)e-^^* 

\6fi^-7fi + 2 ^ 2-3At 



3/i 

Ajt 



(25) 



The cross-correlations between random particle-displace- 
ments perpendicular to qi2 are independent of the Weis- 
senberg number Wi: 



{yim2it))^{smh{t))^\{e 



-Aat 



e-^^*) . (26) 



Both correlation functions are plotted in fig. [6l 

As indicated in fig. [21 the relaxation processes of the 
displacements of the two particles along the same spatial 
direction can be decomposed into a parallel and an anti- 
parallel translation. In the present case the eigenvalue 
Ai (A2) in eqs. (fTO)) corresponds to the parallel relax- 
ation modes along (perpendicular to) the connection vec- 
tor qi2. For the parallel motions, the signs of the particle 
displacements are always equal for both particles (-I-, -I- or 
— , — ), whereas for the anti-parallel ones the correspond- 
ing displacements have opposite signs (-I-, — or — , -I-). As 
a consequence the product of the two displacements is 
always positive for the parallel case and negative for the 
anti-parallel case, as indicated by the prefactors of the 
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FIG. 6. The cross-correlation functions between the two 
beads along parallel directions are shown for <j!> = and 
d = 4a. The minimum of (xi(0)a;2(i)) = {i2(0)ii(t)) in- 
creases with the Weissenberg number Wi and is always deeper 
than the one of the function {yi(Q)y2{t)) = {z\{Q)z2{t)) , which 
is independent of Wi. 



corresponding contributions in the correlation functions 
in eq. ([25)) and eq. (f26|) . As the magnitudes of the anti- 
parallel translations are larger than the parallel ones, the 
superposition of the two different modes is negative and 
has a pronounced minimum at a time t' ss r, as shown 
in fig. [6l The value of t' in case of eq. ([25]) depends only 
weakly on the distance d between the two traps. For the 
cross-correlation (yi (0)^2(0) s-'^ analytical expression for 
t' and its magnitude can be given: 



{vmut')) 




(27) 
(28) 



For finite values of the Weissenberg number Wi 
the cross-correlation of the longitudinal displacements, 
(ii(0)i2(t)), also includes contributions describing re- 
laxation processes perpendicular to qi2 with the eigen- 
values A2 and A4. These additional contributions cause 
larger shear-induced corrections to the cross-correlation 
in eq. ([25l) than for the autocorrelation function in 
eq. ()20ap . which can be seen by comparing the devia- 
tions between the dashed-dotted and the dashed curves 
in fig. [3] and in fig. [G] 

For a single particle in a linear shear flow shear-induced 
cross-correlations between displacements along perpen- 
dicular directions in the shear plane were found, which 
are proportional to the Weissenberg number [3^ . In the 
two-particle system, one also obtains cross-correlations 
between the displacements of two different particles along 
orthogonal directions in the shear plane as described by: 



FIG. 7. Shear-induced cross-correlations between orthogonal 
fluctuations of different particles in the case, where qi2 is par- 
allel to the streamlines, for Wi =1/2 and different distances. 



-A/t 



Wi 



4 V2 + 3At 2-3^4 
Wi 



(y2(0)ii(t)> = — (e-^^* + e-^^*) 

Wi 1^(1 +/i)e-^i* ^ (l-^)e~^3t 
~ 2/1 



2 + 



(29a) 



(29b) 



Both functions are plotted for different values of the trap 
distance in fig. [7] and one can see that the magnitudes of 
the correlations are larger for smaller values of d. Note 
that in general the particle index in the previous equa- 
tions can be interchanged: (.Ti(0)j/2(t)) = {x2{0)yi{t)) . 

The ratio between the two shear-induced cross-corre- 
lations at t = is independent of the Weissenberg num- 
ber Wi and given by: 



-9a 



(5ii(0)yi(0)> 



(30) 



This relation is reasonable for a sufficiently large ratio 
d/a. It might serve in an experiment as a consistency 
check, since the bead radius a and the particle distance 
d are usually well known. If the two beads come close 
to each other, additional efltects related to their shear- 
induced rotation may come into play. This is in general 
a limitation of our Langevin model @ for point-like par- 
ticles with an effective hydrodynamic radius. In Ref. Q 
particle rotations caused by an external torque have been 
taken into account in a Langevin model for two trapped 
particles in a quiescent fluid. This work describes a cou- 
pling between translational and rotational particle mo- 
tions, which was confirmed by experiments |53| . 

The asymmetry in eqs. (|29p with respect to time, 
t ^ —t, has a similar origin as explained above for the 
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shear- induced one-particle correlations ([2T|) . The loca- 
tion of the minimum is again mainly determined by the 
relaxation time r of the beads in the harmonic potentials. 

Note, that for all described correlation functions the 
single particle case is recovered in the limit d — )■ oo, which 
corresponds to /.j — > 0, for example: 



lim(yi(0)ii(t)) = ^ 



1 



2i 

T 



(31) 



All the correlation functions presented in this subsec- 
tion have already been measured in experiments [13, H^l ■ 



B. Perpendicular case: (/> = 7r/2 

For = 7r/2 the time-dcpendcncc of the correlation 
functions and their magnitudes are changed compared to 
the parallel orientation i/i = 0. This subsection focuses 
on these differences. 

The formulas (fT9|) for the eigenvalues of the matrix 
were derived for the identity Rq = Q. This allows the 
determination of exact analytical formulas for the corre- 
lations. In the present case, R7r/2 has to be determined 
numerically as well as the relaxation rates and the ex- 
act functions Cki{t). However, if we use the approxima- 
tion = Q = d/2(0, 1,0,0,-1,0), which is valid for 
small Weissenberg numbers, we obtain analytical formu- 
las for the correlations that also describe quantitatively 
the characteristics of the functions Cki{t) at large values 
of Wi. Within this assumption eqs. (fT9)) remain, but the 
eigenvalues exchange their meanings: Now Ai and A3 be- 
long to the longitudinal fluctuations in the y direction, 
whereas A2 and A4 describe the transversal motions in 
X and z direction, perpendicular to qi2. We also com- 
pare the numerical solutions with direct simulations of 
the Langevin equation ([4]). 



(ii(0)5i(i)) 

■(^i(0)ai(i)> 

■(ii(0)ii(t)> 




FIG. 8. Autocorrelations for </> — n/2, d — Aa and Wi — 
1/2. {xiiO)xiit)), (yi(O)yi(t)) and (2i(0)2i(f)> are differ- 
ent, because of the different relaxation rates and additional 
Wi^-contributions in x direction. One has {xi{0)xi{t)) = 
(2i(0)2i(t)) only for Wi = 0. 



case, because they describe positional fluctuations ei- 
ther parallel or perpendicular to qi2. Only for a van- 
ishing Weissenberg number one obtains (ii(O)xi(i)) = 
(zi(0)zi(<)), but for a finite shear rate the correlation 
functions are all different - for the approximations in 
eqs. (j32[) as well as for the numerical solutions shown in 

fig.El 

The shear-induced cross-correlations between particle 
displacements along orthogonal directions in the shear 
plane behave qualitatively similar as in the case (f) = 0. 
For the approximation R7r/2 = Q they are described by 
the formulas 



1. One-particle correlation 

For the perpendicular configuration the one-particle 
autocorrelations are given for the approximation R^/2 = 
Q by the following expressions: 



(xi(O)xi(i)) = 




-A2t _j_ 








Wi^ 


/(2m 


-l)e 


-A4t 


e-A3t \ 








— 5/.i 


+ 2 + 


2-3/ij 




Wi^ 


/(2/i 


■f l)e- 


-X2t 


e-Ait X 


(32a) 






+ 5m 


+ 2 


2 + 3^7 ' 


(yi(O)yi(t)) = 


4 ^ 


Alt _|_ 


g-A3t) 




(32b) 


(ii(O)zi(t)) = 


4 ^ 


\2t 


g-A.t) 




(32c) 



While the correlations {yi{0)yi{t)) and (zi(O)zi(t)) of the 
displacements perpendicular to the flow lines were iden- 
tical in the case = 0, they are different in the present 



(5i(0)yi(t)) 



Wi 

Wi 



-Alt 



-Agt 



2 + 3m 2-3/i 



(yi(0)5i(i)) = ^(e-^^*-e-^^*) 

Wi f{l + 2^)6-^2* ("1 _ 2M)e-^i* 
^ "2 



in 



2 + 3/i 



2 - 3/1 



(33a) 



(33b) 



The magnitudes of the different contributions in eqs. (j33p 
changed compared to the expressions in eqs. pT|) . which 
influences the positional probability distribution P(r) of 
one particle. The dependence of the inclination angle 
6 on the distance d is weaker in the present case than 
for (/) = 0, see fig. [13] in the next section. Possible con- 
sequences of this difference for the dynamics of beads- 
spring models for polymers in shear flows are discussed 
in the concluding remarks. 

The magnitudes of the correlations functions given by 
eqs. p3p are smaller than those determined numerically 
for the exact distance R7r/2> which are plotted in fig.IH In 
the same figure these exact solutions arc compared with 
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FIG. 9. Comparison of the numerically determined one- 
particle cross-correlations (lines) with the results from the 
direct simulation of eq. Q (circles, squares). The error bars 
are smaller than the symbols. Parameters: d = 4a, Wi = 1/2. 



the results of a Brownian dynamics simulation of eq. (g]), 
where 2*10^ ensemble averages were made. This example 
illustrates the good agreement between both approaches, 
which is not affected by the choice of the Rotne-Prager 
tensor instead of the Oseen tensor in the simulation. The 
analytical expressions ([33)) turn out to be a good approx- 
imation for the parameters d > 8a and Wi < 0.1. Espe- 
cially for smaller values of the trap distance d deviations 
in the relaxation times occur. 



2. Inter-paHicle correlations 

The cross-correlations of the random displacements be- 
tween the two particles in the same direction are given 



for R 



■■it/2 



Q by: 



{Xi{0)x2{t)) = 






g-A4t) 




Wi^ 




2^)6 




-As* 


Wi^ 




— 5/i 
2/i)e 


+ 2 2- 

-A,t 


- 3/Lt 
-Alt 








+ 2 2- 


f 3Ai 






Alt _ 






(il(0)z2(t)) = 




Aat _ 







(34a) 
(34b) 
(34c) 



Similar to the one-particle correlations in sect. IIIIB II 
and in contrast to the case (/) = in sect. IIII A21 all 
three correlation functions are different for finite values 
of the Weissenberg number. This is indicated by the 
expression in eqs. (j34p as well as by the numerical corre- 
lation functions displayed in fig. 1101 and they are again 



FIG. 10. Inter-particle cross-correlations in the orthogonal 
case for d ^ 4a and Wi = 1/2. {xi{0)x2{t)) , {j/i(0)?/2(t)> 
and (21(0)22(4)) are all different and anti-correlated with pro- 
nounced minima. One obtains (xi{0)x2it)) — {zi{0)z2(t)) for 
Wi = 0. 



all anti-corrclated. 

The most striking difference between the case 4> = 
and (j) = 7r/2 is found by looking at the shear-induced 
inter-particle correlations between orthogonal directions 
in the shear plane, which are given for = Q by the 

following expressions: 



(yi(0)i2W) 
Wi 



Wi 

-Wi 



-Alt 



-Aat 



f 3/i 

-Alt 



4/i ^ 
(1 + 2/i)e-^2t 

2 + 3/i 



+ e 



2-3/i 

-A3t\ 



(1 - 2^)e-^-i* 
2-3/1 



(35a) 



(35b) 



A Taylor expansion of the correlation function 
{yi{0)x2{t)) with respect to small values of /x reveals the 
difference between the parallel and the perpendicular ori- 
entations. For (/) = one obtains up to the linear order 

of fi, 



(36) 



(yi(0)i2W> ~ /t (^3 + 4^ + 6^ j e-*/^ , 

whereas for = 7r/2 one obtains, 

(yi(0)i2(t)) ^ /t (^3 + 2^ + 6^^ e-*/^ . (37) 

The first expression has only one cxtremum as a function 
of time. The different prcfactor of te"*/"^ in the second 
expression is the origin of an additional extremum in the 
case (j) = 7r/2. 

For Q = d/2(0, 1, 0, 0,-1, 0) the angle between the vec- 
tor connecting the resulting mean positions of the parti- 
cles, (ri2), and the y axis increases as a function of Wi. In 
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FIG. 11. Shear-induced inter-particle correlations in case of 
a tilted trap for a Weissenberg number Wi =1/2 and for the 
two distances d — Aa and d = 8a. The time-dependence is 
different from the case </> = 0, cf. Fig. [T] 



this configuration the numerical solution for {yi{0)x2{t)) 
does not show a second extremum. However, if Q is tilted 
against the flow direction, like 



Q 



2VT 



Wi^ 



= (-Wi,l,0,Wi,-l,0). 



(38) 



the resulting vector (ri2) becomes nearly parallel to the y 
axis and the second extremum of {yi{0)x2{t)) is obtained 
again, as predicted by the approximation in eq. ([571) . The 
result for the tilted trap is shown in fig. [Tl] for two differ- 
ent values of d. Moreover, the extrema are now stronger 
pronounced than predicted by eq. (jST]) . Comparing this 
with fig. El where the results for = are plotted, it can 
be clearly seen, that the correlation {yi{Q)x2{t)) now has 
an additional local maximum, cf. solid and dash-dotted 
line in fig. [71 and that (ii(0)y2(i)) has a minimum at 
small values of the time t, cf. dashed and dotted lines in 
fig.E 



C. Oblique case: = 7r/4 

If the connection vector qi2 between the two potential 
minima is oblique to the flow direction, further aspects 
for the correlation functions may come into play. In prin- 
ciple, the matrix Cki{t) can be calculated for any angle 
<f>. However, we focus on the special case (/) = 7r/4 as an 
example. This has the advantage that analytical expres- 
sions can be obtained under the assumption = Q. 
These approximate solutions remain rather compact and 
may serve as a guide for the qualitative behavior of the 
correlation functions. 

Since the linear shear flow ^ dictates a preferred di- 
rection in the system, the particle fiuctuations are decom- 



posed into eigenmodes parallel and perpendicular to the 
X axis as described in sect. Ill Bl That's why in the oblique 
configuration the longitudinal displacements consist of a 
superposition of the eigenmodes in the shear plane. If 
(j) = 7r/4, the eigenvalues of the matrix M are given for 
R-7r/4 = Q = dV2/4(l, 1,0, -1,-1,0) by the following 
expressions. 



Ai 

A2 
A3 

A4 

A5 
Ae 



2 + 3/i + ^/i^ - 2/iWi 



2t 



2 + S^i - VM^^AiWi 



2t 



1 + M 



2 - 3^* - + 2MWi 



2t 



2 - i^i + + 2AtWi 



2r 



1-M 



(39a) 

(39b) 

(39c) 

(39d) 

(39c) 
(39f) 



While Ai and A4 describe the parallel and anti-parallel 
relaxation modes in the x direction, A2 and A5 belong to 
the y direction, and A3 and Ag to the z direction. In con- 
trast to the parallel case, the relaxation processes along 
the X and y direction in the shear plane are now affected 
by the shear rate. So the corresponding eigenvalues de- 
pend directly on the Weissenberg number. This was not 
the case in the previous section for R7r/2 = Q- As long as 
Wi 7^ 0, the six relaxation parameters are different and 
can partly become complex numbers, causing oscillatory 
contributions to the functions Cm {t) . 

In order to write down the full expressions for the cor- 
relations in a compact form, we introduce the following 
notions similar to eq. (jlSp : 



(40a) 
(40b) 
(40c) 
(40d) 
(40e) 
(40f) 
(40g) 
(40h) 



91,1 = A* 

92,1 = A* 

54,1 = ^^ 

95.1 = A* 

51.2 = A* 
52,2 = A* 
54,2 = Ai 
95.2 = A* 



Wi(Ai - 


A2) 


2Ai 


7 


Wi(Ai - 


A2) 


2A2 




Wi(A5 - 


A4) 


2A4 


7 


Wi(A5 - 


A4) 


2A5 


7 


Wi(3Ai - 


FA2) 


2Ai 




Wi(3A2 - 


HAi) 


2A2 




Wi(3A4 - 


KAs) 


2A4 




Wi(3A5 - 


HA4) 



2A5 



Analogous to the previous subsections we discuss the 
one-particle correlations first. The autocorrelations of 
a single particle show the same behavior as in the case 



11 



<j> = 7r/2. They are different in distinct directions as 
long as Wi 7^ and decay exponentially in time. For 
R-7r/4 = Q the analytical formulas read 



(ii(0)£i(<)> 



8n 



8^i 



8/i 



8^l 



(yi(o)yi(0> 



8(/i-2Wi) 8(Ai-2Wi) 



(41a) 



(41b) 



8{iJ + 2Wi) 8{^l + 2Wi) ' 
(fi(0)5i(0) = i(e-^^*+e-^«*) . (41c) 

However, in contrast to cqs. ([5^ one can see that the two 
autocorrelation functions for the particle displacements 
in the shear plane, namely {xi{0)xi{t)) and (2/1(0)2/1 (t)), 
include the four corresponding relaxation modes and de- 
pend in a complex way on the Weissenberg number. 

The one-particle cross-correlations between orthog- 
onal positional fluctuations in the shear plane, like 
{xi{0)yi{t)) , are purely shear-induced in the parallel and 
in the perpendicular case. They are linear functions of 
the parameter Wi, which means, that they vanish in the 
limit of zero shear rate. This is different in the orthog- 
onal configuration. For <j> — 7r/4 and R7r/4 = Q the 
correlation functions are given by the expressions: 



(£i(0)yi(0> 



gi.ie 



-Alt 



.92,16 



SVm^ - 2/iWi 8Vm^ - 2/^Wi 



54, le 



55, le 



-A5t 



SVm^ + 2^Wi 8Va*2 + 2/iWi ' 



(42a) 



(yi(O)ii(t)) 



51,26 



-Alt 



52,26 



-Aat 



B^Ai^ - 2^Wi 8 - 2/1 Wi 



54,26 



-Alt 



.95,26 



-Ast 



8v/m^ + 2/iWi 8v//i2 _^ 2/iWi 



(42b) 



In the limit of very large particle distances, when — >■ 0, 
(ii(O)yi(t)) and (?/i(0).Ti(t)) resemble the one-particle 
case given by eq. ([31]l. whereas in the limit Wi — > the 
two functions become equal but do not vanish in contrast 
to the previous subsections. Instead of the approximate 
expressions (|42|) we display in fig. [12] the full numeri- 
cal solution for these correlation functions. If Wi = 0, 
the identity {xi{0)yi{t)) = (j/i(0)ii(i)) is obtained, cf. 
dotted line in fig. [T^ For larger values of Wi the two 
functions resemble the curves shown in fig. |4l 

According to eq. ()22p and cq. (|23|) the inclination 
angle 9 of the elliptical positional probability distribu- 
tion is determined by the static single-particle correla- 
tion functions, which depend on the trap distance d. In 
figfT51 tan(6') is shown as a function of the parameter 
fi — 3a/(4(i) for different trap configurations. The an- 
alytical expression from the parallel case (solid line) is 



X 10 ' 



(ii(0)yi(t)>,Wi = 
(5i(0)yi(t)>,Wi = 0.05 
(:!/i(0)Si(t)),Wi = 0.05 




FIG. 12. Single-particle correlations between orthogonal dis- 
placements in the shear plane for different Weissenberg num- 
bers, = 7r/4 and d = 4a. In contrast to the parallel or 
the perpendicular case, (5;i(0)yi(t)) and {yi{0)xi{t)) remain 
finite for Wi = 0. 



■(t> = 

■ 4> = 7r/2, tilted 
(j) = tt/2 
4) = tt/4 




0.05 



0.1 



0.15 



0.2 



0.25 



FIG. 13. The inclination angle of the single particle distri- 
bution, as given by eq. (|22[) . is shown as a function of the 
parameter jj,. For Wi =1/2 the exact results for the different 
setups described in the text are compared. 



compared with the full numerical solutions for the other 
setups. For (p = tt/2 (circles) the monotonous decrease of 
tan(6') is weaker than for = 0. For the tilted configura- 
tion given by eq. (|38p has a broad maximum (triangles), 
whereas in the oblique case, = 7r/4, the inclination an- 
gle is increasing continuously with increasing fi (squares). 

The inter-particle cross-correlations in parallel direc- 
tions arc again anti-corrclated in the present case. For 
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FIG. 14. Inter-particle cross-correlations in the oblique case 
for 4a and Wi = 1/2. In the limit Wi = 0, {xi{0)x2{t)) 
and (yi(0)j/2(t)) become equal. 



FIG. 15. Cross-correlations in the oblique case for d = and 
different Weissenberg numbers. The amplitude of {yi{0)x2{t)) 
is increasing much stronger with increasing values of Wi than 
the amplitude of (ii(0)j/2(i)). 



by: 

the approximation R7r/4 = Q, wc obtain: (5;i(0)j/2(i)) = 9i,i'^ ^ 92,ie ^ 



8 V/i2 + 2^jWi Sv/a*^ + 2AiWi 
(xi(0)x2(t)) = — + {yi{0)x2it)) ^ 



Sa^ 8/^ BVm^ - 2/^Wi 8 - 2/iWi 

(43a) _ -9^.2^"'^* ^ + g^-2'^''°* ^ . (44b) 



8(/i — 2Wi) 8(/i — 2Wi) illustrated by the numerical results for the correct 

in fig.[T5l it is remarkable that {xi{0)y2{t)) does not 
(43b) change very much with the Weissenberg number, while 



54,26-^** .95,26"^=* 



8{n + 2Wi) 8(/i + 2Wi) ' ^^le amplitude of {yi{0)x2(t)) increases with Wi. Their 

(z (0)z (t)) — - (e^^^* — e^'^"*) (430) asymptotic behavior for Wi — ^ is similar to the single- 

4 ^ ' ' particle correlations given by eqs. (|42|) . 



The eqs. (|43ll indicate that these correlations are different 
from each other, similar to the perpendicular case. The 
correlation functions for the correct Rt^/^ are shown in 

fig- El 

The cross-correlations (ii(0)y2(i)) and {yi{0)x2{t)) 
show the strongest dependence on the orientation of the 
trapped particles with respect to the flow direction. If 
(j) = 7r/4, they have one minimum each, which is different 
from the case = and from the tilted configuration 
psp . Within the approximation, the functions are given 



IV. CONCLUSION 

In this work we investigated the dynamics of two Brow- 
nian particles, each trapped by a harmonic potential and 
exposed to a linear shear flow. The one-particle and the 
inter-particle positional correlation functions, which can 
be measured in an experiment, were calculated by solving 
an appropriate Langevin model. We discussed the corre- 
lations in detail as a function of the distance between the 
two traps, as a function of the Weissenberg number Wi, 
and for three different configurations, where the vector 
connecting the two potential minima was either paral- 
lel, perpendicular, or oblique with respect to the flow 
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direction. This relative orientation strongly affects the 
time-dependence of the correlation functions. For the 
parallel configuration exact analytical expressions were 
presented. Otherwise, we provided numerical solutions 
and analytical approximations for the correlation func- 
tions. 

Although the stochastic forces in our model were as- 
sumed to be uncorrelated along orthogonal directions, 
we found a coupling between perpendicular particle dis- 
placements caused by the shear flow, similar to that in 
Ref. [33^ . The resulting shear- induced cross-correlations 
depend linearly on the Weissenberg number, Wi, and 
occur also between orthogonal fluctuations of different 
particles. These inter-particle correlations have zero, 
one, or two extrema as a function of time, depending 
on the particle configuration. Besides generating new 
cross-correlations, the shear flow causes a contribution 
proportional to Wi^ in the correlation functions of parti- 
cle fluctuations along parallel directions. 

Due to the hydrodynamic interaction between the two 
particles the magnitudes of the one-particle correlations 
are enhanced with decreasing trap distance d, while in 
the limit of lar ge d istances the single-particle results pre- 
sented in Ref. [33| were recovered. Moreover, we found 
a significant impact of the parameter d on the positional 
probability distribution of each particle in the shear flow. 
The shape of the elliptical distribution is tilted and 
stretched when the two particles approach each other in 
the parallel configuration. The same effect is observed 
when the shear rate is increased. So the shear-flow effects 
are enhanced by the presence of a second particle. In the 
oblique configuration the opposite effect is observed. 

The shear-induced cross-correlations investigated in 
this work are of the same origin as the correlations be- 
tween orthogonal fluctuations of a single freely floating 
particle as discussed in Ref. [13, HH- However, if the 
Brownian particles are trapped, a direct experimental de- 



tection of these correlations becomes possible. This has 
been achieved recently in Ref. (33 | , where two polystyrene 
latex spheres were trapped by optical tweezers and ex- 
posed to a linear shear flow in a special microfluidic de- 
vice. With this setup the cross-correlations of the po- 
sitional fluctuations for the parallel case, as shown in 
fig. [71 were measured directly for the first time. More- 
over, we predict an additional cxtrcmum for the shear- 
induced inter-particle correlations, if the two particles are 
trapped perpendicular to the fiow lines, cf. fig. [TT] 

The two hydrodynamically interacting beads are 
treated as point-particles. In forthcoming works, the 
presented results on the shear-induced correlations are 
extended by taking the finite particle extension into ac- 
count. Preliminary investigations show that the effects of 
particle rotations do not change the major trends in this 
work [5^ , and for small Wi the rotation can be neglected 
anyhow. 

The present article provides insight on the Brownian 
motion of two trapped particles that might be useful for 
the analysis of the stochastic dynamics of a bead-spring 
model for polymers too, where the Brownian particles are 
connected along a chain and fiuctuate around some mean 
distance to their next neighbors. The fluctuations along 
and perpendicular to the connection vector between two 
neighboring beads may exhibit similar correlations as for 
the three model configurations investigated in this work, 
cf. fig. [TJ In this spirit, a profound analysis of the 
stochastic motion of a bead-spring model in a linear shear 
fiow is an interesting task and may contribute further to 
the understanding of polymer dynamics. 
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